Approximate thermodynamic structure for driven lattice gases in contact 
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For a class of nonequilibrium systems, called driven lattice gases, we study what happens when two 
systems are kept in contact and allowed to exchange particles with the total number of particles 
conserved. Both for attractive and repulsive nearest-neighbor interactions among particles and for 
a wide range of parameter values, we find that, to a good approximation, one could define an 
intensive thermodynamic variable, like equilibrium chemical potential, which determines the final 
steady state for two initially separated driven lattice gases brought into contact. However, due to 
nontrivial contact dynamics, there are also observable deviations from this simple thermodynamic 
law. To illustrate the role of the contact dynamics, we study a variant of the zero range process and 
discuss how the deviations could be explained by a modified large deviation principle. We identify an 
additional contribution to the large deviation function, which we call the excess chemical potential, 
for the variant of the zero range process as well as the driven lattice gases. The excess chemical 
potential depends on the specifics of the contact dynamics and is in general a-priori unknown. A 
contact dependence implies that even though an intensive variable may equalize, the zeroth law 
could still be violated. 

PACS numbers: 05.70.Ln, 05.20.-y 



I. INTRODUCTION 

Equilibrium systems, whieh satisfy detailed balance 
and therefore do not have any particle or energy cur- 
rent, are based on a well founded thermodynamic the- 
ory. Studies of equilibrium systems start with the zeroth 
law which is the corner-stone of equilibrium thermody- 
namics. The zeroth law states that, in equilibrium, there 
exists a set of intensive variables, each of which being eon- 
jugate to a corresponding extensive conserved quantity, 
and these intensive variables equalize when two systems 
are kept in contact and allowed to exchange the conserved 
quantities. Speeifieally, an equilibrium system in contact 
with a reservoir is characterized by the familiar Boltz- 
mann distribution where the probability of a microstate 
C is given by P{C) - cxp[-/3{iJ(C) - ^iN}] with H{C) 
the internal energy of configuration C, N the number 
of particles in the system, and the intensive variables (3 
and fj. being the inverse temperature and the chemical 
potential of the reservoir, respectively. These variables /3 
and fi are conjugate to energy and particle-number of the 
system. Moreover, in equilibrium, there is a class of gen- 
eral fluctuation-response relations, collectively called the 
Fluctuation Dissipation Theorem (FDT), which relate 
the response of a system upon change of an intensive vari- 
able (e.g., chemical potential) to the fluctuation in the 
corresponding extensive variable (e.g., particle- number). 

One could inquire whether there can be a similar ther- 
modynamic structure for nonequilibrium systems as well. 
Among the vast class of nonequilibrium systems, one 
ubiquitous subclass is that of systems having nonequi- 
librium steady state (NESS) 0. The systems in a NESS 
have time-independent macroscopic properties whieh is 
similar to that of systems in equilibrium. However, un- 
like in equilibrium, they have a steady current and gener- 



ally cannot be characterized by the Boltzmann distribu- 
tions with an a-priori known energy functions. Perhaps 
not surprisingly, even for this conceptually simplest class 
of driven systems with NESS, there is no well founded 
thermodynamic structure. Intensive studies in this di- 
rection to find a suitable framework for description of 
NESS have been undertaken @-Q- In this paper, we ask 
whether a homogeneously driven many-particle system 
can be characterized in terms of an intensive thermody- 
namic variable whieh equalizes for two systems in eon- 
tact. We address this question using a simple class of 
stochastic models called driven lattice gases. 

Although there have been many attempts to define an 
intensive variable for driven systems in various specific 
contexts [7H10| , there was no general formulation in this 
regard until recently when a prescription to define such 
a variable for driven systems was proposed by invoking a 
hypothesis called the asymptotic factorization property 
[llj . This property has been shown to be satisfied for a 
class of driven systems having short-range spatial corre- 
lations such as the zero range process (ZRP). 

The ZRP is one of the simplest example of driven 
interacting many-particle systems whieh do not satisfy 
detailed balance and have nonequilibrium steady states 
[l2|. Previously, the ZRP has been mainly considered as 
model-system for various mass transport processes and 
has been used to study the phenomenon of condensa- 
tion transition in nonequilibrium systems Il3| . Recently, 
it has been demonstrated [1J| that the systems governed 
by the ZRP has a simple thermodynamic structure where 
a suitably defined intensive variable, like equilibrium 
chemical potential, indeed equalizes for two such systems 
in contact. There is also a corresponding fluctuation- 
response relation between the compressibility and the 
fluctuation in the particle-number, which is satisfied ex- 
actly for a system in contact with a particle reservoir. 



The main advantage of studying these simple models like 
the ZRP was that the steady-state probability distribu- 
tions, which have simple factorized forms, can be cal- 
culated exactly and therefore various features of driven 
systems in contact can be studied analytically. Due to 
the simple form of the steady-state distribution, it was 
also possible to analyze the role of dynamics at the con- 
tact between two systems. 

However for driven systems with nontrivial steady- 
state properties, the situation is expected to be far more 
complex. Here we consider a simple model of a driven in- 
teracting many-particle system known as Katz-Lebowitz- 
Spohn (KLS) model ^^. The KLS model, first intro- 
duced to study fast ionic conductors [17|, is a paradigm 
in nonequilibrium statistical mechanics. The model de- 
scribes a stochastic lattice gas of charged particles which 
is homogeneously driven by a constant externally applied 
electric field. Initially, the primary motivation behind 
studying the KLS model was to understand the non- 
trivial spatial structures and phase transitions in generic 
bulk-driven systems with nonequilibrium steady states. 
Since the introduction, the KLS model has been stud- 
ied intensively. By now, the phase diagram in the plane 
of temperature and electric field strength is quite well 
known, mainly from extensive simulations |15l . ll6Lll8l . ll9l | 
as well as the results from mean-field theory [20[ and 
renormalizatio n-g roup analysis of a continuum version of 
the model 19|, |2l| . However there is still no well founded 
thermodynamic theory for these driven interacting many- 
particle systems. 

In this paper, we explore the "equilibration" between 
two driven lattice gases when they are brought into con- 
tact. Recently we studied the KLS model, with repul- 
sive nearest- neighbor interactions among particles |22| . 
which revealed a simple but approximate thermodynamic 
structure. Here we extend our previous studies to the 
systems with attractive interactions as well. Interest- 
ingly, both for attractive as well as for repulsive interac- 
tions and for a wide range of parameter values, we find 
that, to a very good approximation, there is an inten- 
sive thermodynamic variable, like equilibrium chemical 
potential, which determines the final steady state while 
two systems are allowed to exchange particles. Conse- 
quently, the zeroth law of thermodynamics as well as the 
fluctuation-response relation between the compressibil- 
ity and the fluctuations in particle-number are satisfied 
remarkably well in a wide range of parameter values. 

However there are also observable deviations from this 
simple thermodynamic structure, especially at high in- 
teraction strengths and large driving fields. We explain 
these deviations by expressing the asymptotic factoriza- 
tion property, which was initially proposed by Bertin et 
al. in |ll.ll4ll and later discussed by us for driven lattice 
gases in [22|, in a modified form where contributions to 
the large deviation functions due to the contact dynamics 
arc identified. To illustrate the origin of these deviations. 




FIG. 1. A scliematic diagram of two nonequilibrium steady 
states in contact with the contact region Vi and V2. Particles 
are allowed to be exchanged through the contact region Vi and 
V2- In this process the total number of particles N = N1+N2 
are conserved where A^i and N2 are the number of particles 
in system 1 and system 2 respectively. 



we study the nontrivial role of the contact dynamics us- 
ing first a simple variant of the ZRP and later the KLS 
model. We find that, depending on the various parame- 
ter values, the contact dynamics can amount to an excess 
chemical potential across the contact for the variant of 
the ZRP as well as the KLS model. This excess chemi- 
cal potential is generally a-priori unknown, and in some 
sense arbitrary for an arbitrarily chosen contact dynam- 
ics. Therefore, it may not be always possible to assign, to 
the individual systems, an intensive variable which is in- 
dependent of the contact between the two systems, thus 
accounting for the deviations from the zeroth law. 

Here is a brief outline of the paper. In section II we 
describe the model, in section HI we present the numer- 
ical results concerning the zeroth law for the KLS model 
and possible deviations from the law. In section IV, we 
discuss the excess chemical potential for a variant of the 
ZRP as well as a variant of the equilibrium KLS model. 
In section V, we describe how the large deviation princi- 
ple can be written in a modified form which can capture 
the deviations from the zeroth law and then discuss the 
excess chemical potential for the KLS model. At the 
end, we summarize. In Appendix A, we give proof of 
the ansatz for the steady-state probability distribution 
in the case of the ZRP and, in Appendix B, we discuss 
the fluctuation-response relation for the ZRP as well as 
for the KLS model. 



II. MODEL 

We consider stochastic lattice gases of charged par- 
ticles which are driven out of equilibrium by constant 
externally applied electric fields in the bulk and which 
therefore have spatially homogeneous steady states 15 1. 
Particles move on a discrete lattice and jump stochasti- 
cally from one site to any of its nearest-neighbor sites, 
preferably towards the direction of the external driving 
field of magnitude E. Due to hardcore repulsion among 



particles, a lattice site can be occupied by at most one 
particle. In addition, particles may also interact with 
each other through a nearest-neighbor pair-potential of 
interaction strength K. We define the occupation vari- 
able r]{r) at a site r = {r^,, r^} where 77(r) = 0, 1: if a site 
r is occupied r?(r) = 1 otherwise 77(r) = 0. We consider 
two systems of lattice gases where system 1 and system 2 
are confined respectively in volume Vi and volume V2 (see 
Fig. [1]). When two such systems are brought into con- 
tact, they are connected at a finite set of points V( and 
V^; subsets of Vi and V2 respectively (Vl^V^ < Vl,V'2) 
and; while in contact, they can exchange particles with 
each other. The energy function H of the combined two 
systems is given by 

H = K, Y. V{ri)v{r'i) + K2 ^ r;(r2)r7(ri) (1) 



{ri,ri> 



(r2,r^) 



where (*, *) denotes sum over nearest-neighbor sites with 
ri, ri' G Vi and r2, r2' S V2, Ki and K2 strength of inter- 
actions among particles for the respective systems. Note 
that systems may in general have different microscopic 
dynamics depending on the interaction strengths Ki and 
K2. Constant driving fields, Ei and E2 respectively in 
systems 1 and 2, are applied along x-directions, with pe- 
riodic boundary conditions imposed in both x and y di- 
rections. Driven bilayer systems have been considered 
before where particles can jump from one layer to the 
other at any site [23[ . However here we consider the case 
where particles can jump from one system to the other 
through a very small contact region between the systems. 
We choose jump rates of particles such that they satisfy 
the local detailed balance condition [15[ . A pair of near- 
est -neighbor sites, located at r and r', in a configuration 
C arc chosen randomly and an attempt is made to in- 
terchange the occupation variables where the attempted 
final configuration is denoted by C". Let us denote the 
corresponding transition rate from configuration C to C" 
as w{C'\C). For movements of particles inside the same 
system (i.e., particles not jumping from one system to the 
other), a quantity A{E) = H{C') - H{C) - E{r^ - r^), 
which depends on the driving field E, is defined where 
Tx denotes the x-component of the position vector r and 
H{C) is the energy of the configuration C . The transition 
rate is assigned to be 



w(C"|C)=min{l,e-''^(-^'} 



(2) 



where (3 is the inverse temperature of the heat bath. 
When the chosen pair of sites are such that a jump is 
attempted from one system to the other across the con- 
tact, the transition rate is assigned to be 



w(C'|C)=min{l,e-'3'^(")} 



(3) 



where A(0) = H{C') - H{C). Note that there is no 
field along the bonds connecting the two systems. In the 
simulations, we consider two-dimensional systems {V = 



L X L) with periodic boundaries in both directions. We 
put 13 ~ 1 throughout the paper. 

When El = E2 ~ 0, the jump rates satisfy the de- 
tailed balance condition and the configuration C of the 
combined system has the Boltzmann probability distri- 
bution P{C) ~ exp[-F(C)]. For Ei,E2 ^ 0, there is 
a constant current in the steady state. However, unlike 
in equilibrium, the nonequilibrium steady-state probabil- 
ity distribution in general is not given by the Boltzmann 
distribution with an a-priori known energy function and, 
except for a few cases (e.g., with only hard-core inter- 
actions), the steady-state probability distribution is not 
known. When two systems are brought into contact, they 
can exchange particles with the total number of particles 
N = N\^ N2 conserved where A^i = X^r gv '^(''i) ^"^^ 
N2 ~ X^r (=v ^(^"2) are number of particles in system 
1 and system 2, respectively. In the sections below, we 
consider the cases where the conserved quantity is the 
number of particles and therefore attempt to define an 
intensive variable, called chemical potential in analogy 
with equilibrium, which is conjugate to the conserved 
particle- number. 



III. NUMERICAL RESULTS 
A. The Zeroth Law 

First let us describe the zeroth law in the context of 
equilibrium systems. We consider three systems and per- 
form the following three thought-experiments which are 
schematically presented in Fig. [51 In the first experi- 
ment, system 1 and system 3 are brought into contact 
and allowed to exchange particles with each other. Sys- 
tem 1 and system 3 eventually equilibrate and reach a 
final equilibrium state with constant average densities 
rii and n^ respectively. In the second experiment, sys- 
tem 2 and system 3 are separately brought into contact 
and allowed to exchange particles. In this case, the ini- 
tial density of system 2 is tuned such that system 3 has 
the same final density as that of system 3 in the first ex- 
periment. Let us denote the final equilibrium densities 
for system 2 and 3 in the second experiment as 712 and 
71-3 respectively. Now in the third experiment, system 1 
and system 2 with initial densities ni and n2 respectively 
are brought into contact and allowed to exchange parti- 
cles. One could ask what the final densities in this case 
would be. The zeroth law of thermodynamics provides 
the answer that the final densities will not change any 
more and will be exactly equal to the respective initial 
densities. Thus the zeroth law allows us to assign to an 
equilibrium system an intensive thermodynamic variable 
called the chemical potential which equalizes for two sys- 
tems in contact. Note that, in Fig. [2] if one compares 
the density profiles of two systems being in the same col- 
umn, the corresponding average density profiles would be 



exactly same. 
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FIG. 2. Schematic diagram to test the zeroth law of ther- 
modynamics for equilibrium systems using a sequence of 
three thought-experiments where three systems (for simplic- 
ity, shown in two dimensions) are separately kept in contact 
with each other and allowed to equilibrate. Experiment 1: 
System 1 and system 3 in contact. Experiment 2: System 2 
and system 3 in contact. Experiment 3: System 1 and system 
2 in contact. 



We use the above strategy to test the zeroth law for 
systems in nonequilibrium steady states. We perform 
the same set of numerical experiments as described 
above, but now with two of the systems driven out of 
equilibrium due to external driving fields present in the 
bulk of the individual systems. Interestingly, similar to 
the equilibrium case, we observe that, for various values 
of interaction strengths, driving fields and densities, the 
zeroth law is quite well satisfied, i.e, if two driven lattice 
gases are separately "equilibrated" with a common 
system with a fixed density, they will also equilibrate 
amongst themselves. Two such examples are given below. 

1. Example for driven systems with attractive in- 
teractions - In Fig. [3] we consider three systems, a 
driven system 1 with A' = — 1, _B = 6, a driven system 

2 with K = —0.75, E = 4 and an equilibrium system 

3 with K ^ 0, E = Q. First, system 1 with density 
ni ~ 0.60 and system 2 with density n2 ~ 0.29 arc 
separately "equilibrated" with system 3 with a fixed 
density n^ ~ 0.15. We then find that system 1 and 
system 2, with the initial densities ni ~ 0.60 and 
712 — 0.29 respectively, "equilibrate" with each other 
such that, to a very good approximation, the respective 
final steady-state values of densities n[ ~ 0.60 and 
712 — 0-29 remain almost unchanged. 

2. Example for driven systems with repulsive in- 
teractions - In Fig. 21 we consider three systems, a 
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FIG. 3. Numerical experiments to test the zeroth law. Top 
panel: Average densities as a function of position co-ordinates 
rx and ry. Bottom panel: Cross-sections (along 2;-direction) 
of the density profiles. Driven system 1: K — —1, E — 6 (top 
density profiles). Driven system 2: K = —0.75, E — 4 (mid- 
dle density profiles). Equilibrium system 3: K — 0, E — 
(bottom density profiles). All systems considered here are 
with same volume V' = 50 x 50. Experiment 1: System 1 
with density n± ~ 0.60 (top grey profile) "equilibrated" with 
system 3 with density na ~ 0.15 (bottom red profile). Exper- 
iment 2: System 2 with density 722 — 0.29 (middle magenta 
profile) "equilibrated" with system 3 with density na ~ 0.15 
(bottom green profile). Experiment 3: System 1 with den- 
sity n'l (top black profile) "equilibrated" with system 2 with 
density n'2 (middle blue profile) where n'l ~ m and n'2 ~ 712. 



driven system 1 with /i = 3.75, A = 6, a driven system 
2 with K = 1.5, E = 5 and an equilibrium system 3 
with K = 0.75, £■ = 0. System 1 with density m ~ 0.40 
and system 2 with density 712 — 0.57 are separately 
"equilibrated" with system 3 with a fixed density 
713 ~ 0.76. In Fig. IH one could see that system 1 and 
system 2 with the above initial densities ni ~ 0.40 and 
7i2 — 0.57 respectively equilibrate with each other with 
the respective final densities 77'^^ ~ 0.39 and 7^2 — 0.58. 
Also in this case, the zeroth law is satisfied to a good 
approximation where final densities n'l and 712 remain 
almost same as the initial densities. 

Interestingly, as seen in Figs. [3] and [H two systems 
in contact have homogeneous density profiles even when 
one or both of them may be driven. The driven systems 
arc indeed far away from equilibrium since the numeri- 
cal values of currents in the nonequilibrium steady states 
considered in Figs. [3] and [3] are near to the respective 
maximum values of currents (data not shown). More- 
over in the top panel of Fig. |4l system 1 with density 
711 — 0.40 has a homogeneous disordered state in con- 
trast to an ordered state for the corresponding equilib- 
rium system. The equilibrium system, with the same 
interaction strength K — 3.75 and the same density, has 
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FIG. 4. Numerical experiments to test the zeroth law. Top 
panel: Average densities as a function of position co-ordinates 
rx and Vy. Bottom panel: Cross-sections (along 3;-direction) 
of the density profiles. Driven system 1: K = 3.75, E = Q 
and 1/ = 32 X 32 (bottom density profiles). Driven system 2: 
K — 1.5, E = h and V = 40 x 40 (middle density profiles). 
Equilibrium system 3: iC = 0.75, E = Q and F = 50 x 50 (top 
density profiles). Experiment I: System f with density ni ~ 
0.40 (bottom red profile) "equilibrated" with system 3 with 
density 713 ~ 0.76 (top grey profile). Experiment 2: System 2 
with density 712 — 0.57 (middle black profile) "equilibrated" 
with system 3 with density 713 ~ 0.76 (top green profile). 
Experiment 3: System 1 with density n'l (bottom blue profile) 
"equilibrated" with system 2 with density n'2 (middle magenta 
profile) where n'l ~ ni and n'2 — n2- 
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a symmetry-broken ordered phase with a checkerboard- 
Hke pattern where sub-lattice densities are different [l6[ . 
Importantly, the macroscopic properties like densities 
do indeed depend on the driving field when a driven sys- 
tem 1 is kept in contact with an equilibrium system 2 
with a fixed density 7i2 • This can be seen in the behavior 
of density ni of a driven system as a function of driving 
field El . We consider the driven systems which were pre- 
viously used to test the zeroth law in Figs. [3]and|4l In the 
top panel of Fig. [51 densities of two driven systems with 
attractive interactions, with Ki = —1 and Ki = —0.75, 
have been plotted as a function of the driving fields Ei in 
the respective systems where both the systems arc sep- 
arately kept in contact with an equilibrium system with 
K2 = and a fixed density 712 — 0.15. In the bottom 
panel of Fig. [5l densities of two other driven systems with 
repulsive interactions, with Ki = 3.75 and Ki ~ 1.5, 
have been plotted as a function of the driving fields Ei in 
the respective systems where both the systems are now 
separately kept in contact with an equilibrium system 
with K2 = 0.75 and a fixed density n2 — 0.76. One could 
see that the densities vary quite significantly, almost by 
an amount of 10% or more from the respective equilib- 
rium values depending on the interaction strengths, when 



FIG. 5. Density rii of driven system 1 as a function of driv- 
ing field El when the system 1 is in contact with an equilib- 
rium system 2 with a fixed density 712. Top panel: Densities 
of two driven systems with interaction strengths Ki = — 1 
and Ki — —0.75 plotted as a function of respective driving 
fields El when the systems are separately in contact with 
an equilibrium system with K2 = and 722 — 0.15. Bot- 
tom panel: Densities of two driven systems with interaction 
strengths Ki = 3.75 and Ki = 1.5 plotted as a function of 
respective driving fields Ei when the systems are separately 
in contact with an equilibrium system with K2 = 0.75 and 
712 ^ 0.76. 



the driving field varies from Ei = to a large value 
El >Ki. 

Provided that the zeroth law is satisfied, one can define 
a chemical potential even for a driven system as follow- 
ing. A driven system is kept in contact with an equilib- 
rium one and allowed to reach a steady state. Then, in 
the steady state, one can simply assign the chemical po- 
tential of the equilibrium system to the driven one. For 
nonzero interaction strengths K ^ 0, even the equilib- 
rium chemical potential fj, cannot be calculated directly 
as one does not know the explicit form of fJ.(n) as a func- 
tion of density n. However, for equilibrium system with 
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FIG. 6. Density n vs. chemical potential fi is plotted for a 
driven system with Ki = — 1, _Ei = 6 which is separately kept 
in contact with three equilibrium systems with interaction 
strengths A'2 ~ 0, K2 ~ —0.25, and K2 = —0.5. The collapse 
of the curves agrees well with the zeroth law. 



noninteractiiig hardcore particles K ^ 0, the chemical 
potential /i can be expressed as a function of density n 
where 



M 



In 



(4) 



by using the relation /i = —{ds/dn) where s ~ — [nlnn + 
(1 — n) ln(l — n)] is the equilibrium entropy per lattice 
site. Since the zeroth law is exactly satisfied in equilib- 
rium, the chemical potential for any equilibrium system 
with K y^ can be measured by keeping the system in 
contact with an equilibrium system with K = and then 
assigning a chemical potential for a system with A' ~ 
to that with K ^ 0. In Fig. El we plot the density n 
as a function of chemical potential /i when a driven sys- 
tem with Ki = — 1, i?i = 6 is separately kept in contact 
with three equilibrium systems with various interaction 
strengths A'2 = 0, K2 = -0.25, and K2 = -0.5. Given 
that the zeroth law is satisfied, the various curves for 
density as a function of chemical potential should fall on 
each other. The collapse of curves in Fig. IH] indeed agrees 
quite well with the zeroth law. 

We have studied the dependence of densities on the 
system sizes as well. We find that the densities of two 
systems in contact with each other, when one or both 
systems may be driven, are almost independent of system 
sizes. In Fig. [7] we plot the density n as a function 
of chemical potential fj. for various system sizes, and for 
various values of interaction strengths and driving fields. 
The densities seem to only depend on the interaction 
strengths and driving fields of the respective systems. 

In Figs. [3] and |31 we have tested the zeroth law only 
for a particular set of values of densities. We now study 




FIG. 7. Density n as a function of chemical potential /i for 
driven systems (fixed Ki, Ei) with different values of system 
size Vi in contact with an equilibrium system of noninteract- 
ing hardcore particles (A'2 — 0,E2 — 0) with different values 
of system size V2. Circles correspond to the case when the 
system 1 with Vi = 32 x 32 is in contact with system 2 with 
V2 = 32 X 32. Squares correspond to the case when the sys- 
tem 1 with Vi = 20 X 20 is in contact with system 2 with 
V2 = 100 X 100. For the curves from left to right: (1) two 
curves (red) for systems with A'l — —1, Ei = 6, (2) two 
curves (blue) for systems with Ai — —0.75, Ei = 4, (3) two 
curves (sky blue) for systems with Ai = 1, Ei = 6, and (4) 
two curves (black) for systems with Ai = 2, Ei = 2. 




FIG. 8. Final densities n'l and n2 for two driven systems, 
respectively with K = -1, E = 6 and A = -0.75, E = 4, 
being in contact with each other are compared with their 
respective initial densities ni and 722. The initial densities 
correspond to those obtained by keeping the driven systems 
separately in contact with an equilibrium system of a fixed 
chemical potential /i. If the zeroth law is satisfied, the initial 
and final densities should be exactly the same for a given fi. 



the zeroth law for various other densities where the in- 
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FIG. 9. Final densities n^ and n'2 for two driven systems, 
respectively with K = 3.75, -E = 6 and K = 1.5, E = 5, 
being in contact with each other are compared with their 
respective initial densities ni and 712. The initial densities 
correspond to those obtained by keeping the driven systems 
separately in contact with an equilibrium system of a fixed 
chemical potential /i. If the zeroth law is satisfied, the initial 
and final densities should be exactly the same for a given fi. 



teraction strengths and driving fields are kept fixed. We 
consider driven systems 1 and 2 vifith respective densities 
ni and 712 where the systems arc separately equilibrated 
with an equilibrium system 3 with density n^ and corre- 
sponding chemical potential fi. Then system 1 and sys- 
tem 2 with the respective initial densities ni and 712 are 
equilibrated with each other where they eventually reach 
final densities n'l and n'2. Clearly, if the zeroth law is sat- 
isfied, then the corresponding initial and final densities 
would be exactly the same, i.e., ni = n'l and ^2 = n'j. 
We consider the systems as previously discussed in Figs. 



[3] and 131 In Fig. [SJ we plot ni and n'l for system 1 
with K ~ —1, E = 6 and 712 and n'j for system 2 with 
K = —0.75, £' = 4 as a function of chemical potential ji 
of the equilibrium system 3, with K = E = 0. Similarly, 
in Fig. ini we plot 77,1 and n[ for system 1 with K = 3.75, 
E = 6 and 772 and n'2 for system 2 with K = 1.5, E ^ 5 
as a function of chemical potential fi of the system 3, 
with K = 0.75, E ~ 0. Although, the zeroth law is sat- 
isfied to a good approximation, there are indeed small 
but observable deviations from the law, i.e., up to 5% 
deviations in the final densities from the corresponding 
initial density values. In the next section we discuss the 
deviations from the zeroth law in more detail. 



B. Deviations from the Zeroth Law 

In the previous section, we have seen that, for a large 
parameter range, the driven lattice gases have an approx- 
imate but remarkably simple thermodynamic structure. 
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FIG. 10. In top panel, the cross-sections (along s-direction) 
of the density profiles are plotted. Systems 1, 2 and 3 (all with 
same size V = 50 x 50) have interaction strengths K = —1 
(top density profiles), K — —0.75 (middle density profiles) 
and K = (bottom density profiles) respectively. First, sys- 
tem 1 with density ni and system 2 with density 722 are sep- 
arately equilibrated with system 3 with density 713. Then 
system 1 and 2, with initial densities ni and n2 respectively, 
are equilibrated with each other where the final steady-state 
densities are n'l and n'2 respectively. Top panel (left): For 
equilibrium systems where ni = n'l and 712 = n'2. Top panel 
(right): For the case when system 2 driven with field E2 = 10 
with ni 7^ 71 '1 and 712 "^ n'2. Bottom panel: Difference in den- 
sity 5n2 =712 — 712 is plotted versus driving field E2 for system 
2. The blue line, which is a fitting function (a — be~'^^^) with 
a = 0.036, b — 0.036 and k — 0.78, is a guide to the eye. 



However, it should be noted that there are also larger ob- 
servable deviations from the simple thermodynamic law 
as discussed next. 

In Fig. IIOI we perform numerical experiments similar 
to that discussed before (e.g., see Figs. [3]and|4]) where 
system 1 with density 711 and system 2 with density 712 
are separately equilibrated with system 3 with a fixed 
density n^. Then system 1 and system 2, with the re- 
spective initial densities ni and 712, are equilibrated with 
each other where the final steady-state densities are n'l 



and Ti'2 respectively. In the top-left panel of Fig. [TOl 
six density profiles along x-direction are plotted for equi- 
librium systems 1, 2, and 3. Expectedly, for equilib- 
rium systems, the zeroth law is exactly satisfied where 



ni 



0.79 (top density profiles) and 712 



0.41 



(middle density profiles), i.e., initial densities are equal 
to the respective final densities. In the top-right panel of 
Fig. [ini six density profiles along a;-direction are plotted 
when systems 1 and 3 are equilibrium systems but sys- 
tem 2 is driven with a large field E2 = 10. In this case, 
the zeroth law is observed to be violated significantly 
where ni ~ 0.79 ^ n[ ~ 0.75 (top density profiles) and 
7T.2 ~ 0.34 7^ 712 — 0.38 (middle density profiles), i.e., final 
densities change appreciably as compared to the respec- 
tive initial densities. In the bottom panel of Fig. [101 the 



density difference Sn2 



712 is plotted as a function 



of the driving field £'2 in system 2. 
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FIG. 11. Density ti as a function of chemical potential fi for a 
driven system with Ki = —1, Ei — 6 which is separately kept 
in contact with three equilibrium systems with interaction 
strengths K2 = 0, A'2 = -0.8, and K2 = -0.9. 



One can also observe the deviations from the zeroth 
law which is now studied in a slightly different way as 
follows. We first try to assign a chemical potential to a 
driven system by keeping the system separately in con- 
tact with various equilibrium systems and then compare 
density vs. chemical potential curves. In Fig. Ill[ we 
plot density 7i as a function of chemical potential fi for 
a driven system, with Ki = —1, Ei = 6, which is 
separately kept in contact with three equilibrium sys- 
tems with the following interaction strengths: K2 = 0, 
K2 = —0.8 and K2 = —0.9. If the zeroth law was satis- 
fied, all the curves should fall on each other. However, we 
observe that there are significant deviations which occur 
especially around density n = 1/2. In Fig. [T^ we plot 
n vs. 11 for another driven system, with repulsive inter- 
action Ki = 3.75, El = 6, being separately in contact 
with equilibrium systems with K2 = 0.75 and K2 = 1. 
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FIG. 12. Density 71 as a function of chemical potential ^ for 
a driven system with Ki = 3.75, -Ei = 6 which is separately 
kept in contact with two equilibrium systems with interaction 
strengths K2 = 0.75 and K2 ~ 1. 



We again see deviations from the zeroth law when the 
density approaches n = 1/2. 



IV. ROLE OF CONTACT DYNAMICS: 
RESULTS 



EXACT 



A. Excess Chemical Potential in the ZRP 

We now study a simple class of driven systems, known 
as the zero range process (ZRP) [1^, to understand the 
role of contact dynamics in the context of equalization of 
thermodynamic variables for driven systems. The ZRP 
is defined on a discrete lattice where, unlike the KLS 
model, there is no restriction on the occupation number 
of a site, i.e., the lattice sites can be occupied by more 
than one particle. The jump rate of a particle out of any 
site is assumed to depend on the number of particles on 
the site. 

For simplicity, we consider two one-dimensional rings, 
ring 1 and ring 2, consisting of Li and L2 sites respec- 
tively. The rings are kept in contact with each other 
so that they can exchange particles through the con- 
tact area. In the ath ring {a = 1,2), any site ia 
(ia = 1,2, ...,^0) is occupied with rii^ particles. Any 
configuration C can be specified, by using the occupation 
numbers at all sites of two rings, C = {{rii^} , {ui^}) . Two 
rings are connected at two sites, say ii = 1 and 72 = 1, 
through which particles can be exchanged between the 
rings. The dynamics is defined as follows: a particle at 
site ia can jump only in one direction, say to the nearest 
neighbor site 7 q, -I- 1 in the clockwise direction (therefore 



violating detailed balance), with rate 



UairiiJ = Va 



fairii, 



fai^ij 



with a = 1,2 



(5) 



where fa («»„ ) is a function of occupation number rii^ and 
the factor Va is independent of Ui^ . The form of the jump 
rates is the same as given in Eq. [5] irrespective of jumps 
in the bulk or jumps from one ring to the other. But the 
factors Va can in general be different, i.e., Vq = Wq when 

(c) 

a particle jumps in the bulk of the ath ring and Va = Va 
when a particle jumps at the contact from ath ring to 
the other ring. For the zero range process, the properties 
of nonequilibrium steady states while two systems are 
in contact has been previously studied, but only for a 



special case when Va 



(b) 



Va and vl 



,('=) 



^ 



Here 



we consider the more general case when Va 7^ Va as well 



,(^) 



„(=) 



as u} 7^ U2 : i.e., the factor Va taking different values 
at the bulk and the contact of the two rings. However, it 
would be easy to see that the factor Va taking different 
values in two different rings does not change the steady- 



(b) 



state properties. Therefore, now onwards, we put v\ 
1. Also, for notational simplicity, we denote Vi 



(b) 



and V2 



.,(c) 



For completeness, let us first discuss 



the special case when vi = V2 [14|. The steady state 
probability distribution can be written in a factorized 
form 






L2 



n /i(^n) n /2(^ 



x6iNi+N2-N), (6) 

where A'^i = ^^ ^^^ tii-^ and N2 — X)i Li "'i2 ^^^ the num- 
ber of particles in ring 1 and ring 2 respectively, Zjv is 
the normalization constant. The delta function in the 
above equation takes into account that the total num- 
ber of particles N = Ni + N2 is conserved. Clearly in 
this case, the joint probability distribution P{Ni,N2) of 
particle numbers iVi and N2 can be written in a product 
form 



P{Ni,N2) 



Zi{Ni)Z2iN2) 



•^N 



(7) 



SlI^^^^t=l /"KJ'^Ct" 



where Za{Na) 

Na). As discussed in [lllllJl, in this case, one can define 
an intensive variable /ia = —dlnZa/dNa which equalizes 
when two rings are kept in contact, i.e., /ii = /i2- 

Now we discuss the general case when vi ^ V2 which 
leads to our main points. Interestingly, as shown in Ap- 
pendix A, the steady state probability distribution can 
still be written in a factorized form 






L2 



n/iK)n/2(^ 



xeAiJVieA2JV2j(^^ 



+ N2-N) (8) 



where jli = ln(l/t!i) and /i2 — ln(l/u2), which we call 
excess chemical potentials. Now the joint probability dis- 
tribution P{Ni,N2) of particle-numbers A^i and N2 can 
be expressed as 

Zn 

The macrostate is obtained by maximizing lnP(7Vi, A''2), 
i.e., dlnP{Ni,N2)/dNi = 0. Therefore it straightfor- 
wardly follows that 



dlnZi 
dNi 



A*i = 



dlnZ2 
dN2 



("2 



(10) 



or, in other words, there indeed exists new intensive 
variable n'a = {—d\nZa/dNa +lnva), with a = 1,2, 
which takes the same values for two rings in contact, i.e., 
/^i ~ /^2- -f"'^'" ^^^ special case when vi = V2, the excess 
chemical potentials are equal (i.e., fli — JI2) and drop out 
of Eq. [TO] which then implies that the old variables ^1 
and fj,2 equalize. But in the case when vi 7^ W2, the new 
variable ^^ takes the role of chemical potential which 
then equalize upon contact. This identification of fx'a as 
an intensive variable which equalizes for two systems in 
contact was missed in Ulll^ where it was concluded that 
there is no such equalization when vi ^ V2. Moreover, 
in contrary to the suggestion in [ij], the detailed bal- 
ance condition is still satisfied at the contact even when 

Vl ^ 1)2. 

For a class of systems specified by a particular set of 
parameters {fa{n),Va}, it can be immediately checked 
that the zeroth law is indeed satisfied. In Fig. [TOl we 
perform various numerical experiments similar to those 
discussed in Fig. [3] and S) We choose /a(n) ~ n^"'^ 
with various sets of parameter values for {Sa,Va}- In 
the first experiment, ring 1 with di = 3, vi = 1 and 
density ni ~ 2.3 equilibrated with ring 3 with S2 = 5, 
V2 = 0.5 and density 713 ~ 7.7. In the second experiment, 
ring 2 with Si = 4, wi = 0.75 and density 712 — 3.8 
is equilibrated with ring 3 with 82 = 5, V2 = 0.5 and 
density 71,3 ~ 7.7. In the third experiment, ring 1 with 
61 ~ 3, Vl ~ I and density n'^ ~ 2.3 is equilibrated with 
ring 2 with 62 = 4, V2 = 0.75 and density n'2 ~ 3.8 
where n'l ~ rii and n'2 — n2- In this case one can see 
that the zeroth law is satisfied (see the left panel of Fig. 
[T3)) . We have also checked that values of the densities 
are such that chemical potential ^^ takes equal values 
for two rings in contact. In the right panel of Fig. [TOl we 
again perform three similar numerical experiments but 
now, in the second experiment, the values of the factors 
Vl = 0.85 and V2 = 0.4 are slightly perturbed by an 
arbitrary amount from the earlier values of the factors 
Vl = 0.75 and V2 = 0.5 which were chosen in the second 
experiment of the first set of numerical experiments (i.e., 
corresponding to the left panel of Fig. [T5|) . One can see 
that, in this case, the zeroth law is violated even though 
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FIG. 13. Density profiles for rings 1, 2 and 3 (all with same 
size L = 100). Left panel: Three numerical experiments to 
test the zeroth law in the case of the ZRP. Experiment 1: Ring 

1 ((5i = 3, wi = 1) in contact with ring 3 (52 = 5, «2 = 0.5). 
Experiment 2: Ring 2 (5i = A, v\ = 0.75) in contact with 
ring 3 {82 = 5, 1)2 = 0.5). Experiment 3: Ring 1 [Si — 3, 
1)1 — 1) in contact with ring 2 (82 = 4, U2 = 0.75). The zeroth 
law is satisfied in this case. Right panel: Three numerical 
experiments to illustrate the role of contact dynamics for the 
deviations from the zeroth law. Experiment 1: Ring 1 (5i = 3, 
vi = 1) in contact with ring 3 {S2 = 5, U2 = 0.5). Experiment 

2 with slightly perturbed jump rates at the contact: Ring 2 
(Si = 4, 1)1 = 0.85) in contact with ring 3 {82 = 5, 1)2 = 0.4). 
Experiment 3: Ring 1 {Si = 3, vi = 1) in contact with ring 
2 {52 = 4, W2 = 0.75). The zeroth law is not satisfied in this 
case. 



the factorization property as given in Eq. |8] exactly holds 
in each of the three experiments. This illustrates the 
role of the contact dynamics for the zeroth law which is 
satisfied only for a precise set of jump rates at the contact 
with no arbitrariness is allowed in these rates. 



B. Excess Chemical Potential in the Equilibrium 
KLS Model 

The above situation in the case of the zero range pro- 
cess (ZRP) is very similar to the special case of the equi- 
librium KLS model when one chooses the transition rates 
as given in Eqs. [2] and [3] with no driving in the bulk of 
the individual systems, i.e., Ei = E2 = 0, and a slightly 
modified transition rates at the contact as follows. In 
general, the transition rate for a particle jumping from 
system 1 to system 2 is given by 



w{C -^ C) = ui X min{l. 



=-A(o) 



} 



(11) 



and the reverse transition rate for the particle jumping 
from system 2 to system 1 is given by 



w{C' ^ C) = wa X min{l,e~^("'} 



(12) 



where we put /3 = 1. For ui = U2 = 1, the jump rates 
at the contact are same as given in Eq. |3l Clearly the 
modified transition rates at the contact amount to an 
additional field £'12 = ln(wi/w2), say from system 1 to 
system 2, along the bonds connecting the two systems. 
Note that the field E12 is not a driving field and just 
introduces an extra overall shift in the chemical poten- 
tial of the system 1. The field E12 accordingly modifies 
the energy function from H to H' by introducing an extra 
chemical potential in Eq. [l] i.e., H' = H+E12N1. In this 
case, the detailed balance is satisfied with respect to the 
Boltzmann distribution with the modified energy func- 
tion H'. Consequently, one could effectively think of an 
excess contribution to the equilibrium free energy of the 
system 1 due to the shift in the chemical potential of the 
system 1 by an amount E12 . Now the condition of mini- 
mization of the total free energy F = Fi -I- _F'2 , as given in 
Eq. [To] for the zero range process, would imply that the 
new intensive variables /i^ = {—dlnZa/dNa + lnVa), not 
the variables Ha = —din Za/dNa, equalize where free en- 
ergies of the respective systems are Fi = —lnZi+ jliNi 
and F2 = — In Z2 + /i2-^2 with partition function Z^ = 
'^(j e^^°'^'~'' , the energy function for individual system 
Ha = KaJ2{r^.r'j 'n{^a)v{^'a) (^cc Eq. [J) for a = 1,2. 
Consequently, the zeroth law is satisfied in the case when 
each system is assigned a particular set of values of the 
factors Va- However, if the factors Va are chosen arbi- 
trarily in any of the numerical experiments as discussed 
in the right panel of Fig. [13] in the case of the ZRP, the 
zeroth law would not hold even for the equilibrium KLS 
model with the modified jump rates at the contact. 



ANALYTICAL APPROACHES FOR THE 
DRIVEN KLS MODEL 



In this section, first we discuss how the zeroth law, 
which has been seen to be satisfied to a very good ap- 
proximation in the simulations of the driven KLS model 
for various parameter values, could be explained in terms 
of large deviation principle. Then we discuss how the 
deviations from the zeroth law, also observed in the sim- 
ulations, could be explained by modifying this large de- 
viation principle. 



A. Large Deviation Principle 

In equilibrium, the zeroth law can be derived from 
variational principles, e.g., by maximization of entropy 
of an isolated system or by minimization of free energy 
in the case when the system is in contact with a reser- 
voir. For some nonequilibrium systems, there may be 
a similar principle, called the large deviation principle 
[3, ]2M ■ Then, the zeroth law can be derived from the 
large deviation principle along the same line as in equi- 
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librium in the following way. Let us consider the follow- 
ing scenario. Two systems are kept in contact with each 
other with a particular dynamics specified at the contact. 
The two systems exchange according to the contact dy- 
namics some conserved quantity, say number of particles, 
such that A^i + A^2 = ^ = constant with A^i , N2 being 
the number of particles in systems 1 and 2 respectively 
(schematically shown in Fig. [1]). The quantities Ni, N2 
are considered to be extensive, i.e., proportional to the 
volume Vi , V2 of system 1 and 2 respectively. Now we are 
interested in large deviations of A'^i and N2 and assume 
that the probability of large deviation P{Ni, N2) in the 
quantity A^i and N2 is given by 



P(iVi,7V2) 



'Vi/i(ni)g-y2/2(n2) 



(13) 



in the limit oi Ni,N2,Vi,V2 ^ 1 keeping respective den- 
sities ni = Ni/Vi and n2 = N2/V2 finite with the nor- 
malization constant exp[— F(iV)]. The above equation is 



the statement of the large deviation principle [2, l2j| and 
the functions /i(ni), f2{n2) are called the large devia- 
tion functions for the corresponding systems. The sign 
'~' implies equality in terms of logarithm of the respec- 
tive quantities and Eq[T3]can be written more rigorously 
as 

- \nP{Ni,N2) = Vi/i(ni) + 1/2/2(^2) - F{N) 

+e{Ni,N2){U) 

where e(A^i,iV2)/lnP(iVi, A^2) -> as A^i,iV2 -^ 00. 
Note that, in writing Eqs. [13] and [141 '^e have assumed 
that the correlation between system 1 and 2 can be ne- 
glected as a boundary-effect in the limit of large volume. 
The macroscopic stationary state {Ni,N2}, under the 
constraint A^i + N2 ~ constant, can be determined by 
maximizing lnP(A^i,iV2), i-c, d\nP{Ni,N2)/dNi = 
which gives 



dnj. \dn2). 



(15) 



where n\ ~ Nl/Vi, rij = N2/V2 and /i is the chemi- 
cal potential which takes same values in the steady state 
when system 1 and 2 are kept in contact. Note that 
the validity of Eq. [151 which implies the existence of an 
intensive thermodynamic variable, follows from the as- 
sumption of a large deviation principle of Eq. 1131 The 
above arguments are quite general and could be valid 
irrespective of a specific nature of the dynamics one con- 
siders in a particular problem. Therefore Eq. [151 can be 
equally applicable to an equilibrium as well as a nonequi- 
librium steady state. Another consequence of equaliza- 
tion of intensive variables is the fluctuation-response re- 
lation between the compressibility and the fluctuations 
in particle- number iVi of system 1, i.e., 



d{N^) 



(n!) {N,)' 



(16) 



when system 1 is kept in contact with a very large system 
2, with chemical potential fi, which can be considered as 
a particle reservoir (i.e., iVi <^ N2 and Vi <C V2). For 
more details regarding the above fluctuation relation in 
the context of the ZRP as well as for the KLS model, see 
Appendix B. 



B. Modified Large Deviation Principle 

The exact results for the ZRP and the equilibrium KLS 
model give us an insight into the role of the contact dy- 
namics which can effectively introduce an excess chemi- 
cal potential across the contact region. Consequently, the 
intensive variable for a system which equalizes upon con- 
tact can have different functional forms depending on the 
specifics of the dynamics at the contact and therefore on 
the other system in contact. To describe this situation in 
a more general context where there may be equalization 
of an intensive variable but still the zeroth law docs not 
strictly hold, we write the large deviation probabilities as 
given in Eq. [131 in a modified asymptotic form 

g-Vi/i(ni) -V2/2(n2)g/ilJVi A2JV2 
^(^1,^^2) ^ZTW (^^) 

where /i , /2 are the large deviation functions for a puta- 
tive contact dynamics for which the zeroth law is satis- 
fied, the two additional quantities fli, JI2 can be thought 
of as excess chemical potentials arising solely due to the 
actual contact dynamics which is under consideration 
and for which the zeroth law is not satisfied. Clearly, 
for an arbitrarily chosen contact dynamics, the old in- 
tensive variables /xi = dfi/dni and /i2 = df2/dn2 do 
not equalize. Note that, unlike in the case of the ZRP or 
the equilibrium KLS model, the potentials /ii and p,2 is in 
general a-priori unknown for a given contact dynamics. 
Moreover, for an arbitrarily chosen contact dynamics, it 
may not always be possible to assign, to the individual 
systems, an intensive variable which is independent of 
the specifics at the contact between two systems. In this 
case, even though there can be equalization of some in- 
tensive variable when two systems are kept in contact, 
the zeroth law may not be satisfied as demonstrated in 
the case of the ZRP (see right panel of Fig. [13]) . 

The above modified large deviation principle in Eq. [T7] 
still assumes that two systems in contact have an asymp- 
totic factorization property in the sense that the corre- 
lations between the systems can be ignored in the large 
volume limit. In a special case (e.g., the ZRP or the 
equilibrium KLS model), the chemical potentials jli and 
JI2 can be constant over an entire density range. How- 
ever, in general, the quantities /ii(ni) and jl2{n2) can be 
functions of densities ni and 712. Note that, although the 
large deviation principles in Eqs. [TS] and [17] appear to 
have different forms, they are essentially the same. The 
only difference is that, in the modified form of Eq. I17[ we 
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identify the contribution to the large deviation function 
due to the contact dynamics and separate this contribu- 
tion from that arising due to the bulk of the individual 
systems. 
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FIG. 14. Correlations between various neighboring sites as 
a function of distance when a driven system 1 [K = —0.75, 
E = A, V = bQ X 50) and an equilibrium system 2 {K — 
-0.75, S = 0, V = 50 X 50) in contact with each other. The 
function q? (r^.a) denotes the scaled correlation (see Eq. llSp 
along ath direction between two points located in system i 
and system j where Vj^a being the ath component of relative 
position vector rj between the two points with i,j = 1,2 
and a — x,y (for two dimensional system), Ar7(r) — r;(r) — 
{rj{r)) the fluctuation in occupation variable rj{r). For i ^ j, 
c\" {rj^a = 0) is the correlation between two nearest-neighbor 
sites located across the contact. 



The asymptotic factorization can indeed be a very 
good approximation even in the case of a more compli- 
cated driven lattice gas like the KLS model. To show 
this, we now study various spatial density-correlations 
q" (rj.c) between two points located in the individual 
systems as well as two points located in two different 
systems across the contact. We define the function 



c!f(^,,o) = 



({A7?(rf)}{Ar;(rJ + rj)}) 
V({A7Krf)}2)y^({Ar;(rJ + rj)}2)' 



(18) 



which denotes the scaled correlation along ath direction 
between two points in system i and system j where rj^a 
the ath component of position vector rj with i,j ~ 1,2 
and a = x,y (for a two-dimensional system), A?/(r) = 
77(r) — (?7(r)) the fluctuation in occupation variable ri{r), 
r|^ the position vector of the contact site in system i. 

For i = j, c\f {vi^a) denotes density-correlations, along 
a-axis, between two points both located in the same sys- 
tem i and, for i ^ j, c]j{rj^a) denotes cross-correlation, 
along a-axis, between two points one located in the sys- 
tem i and the other located in system j. In Fig. I14[ 



we have plotted the scaled correlation q" (^j^a) as a 
function of rj^a for a driven system with Ki = —0.75, 
El ~ A and an equilibrium system with K2 = —0.75 
where the systems are kept in contact, both at den- 
sity n = 1/2. One can see that the amplitude of the 
cross-correlations among nearest-neighbor sites located 
in two different systems across the contact are very small, 
almost an order of magnitude smaller as compared to 
those among neighboring sites in the individual systems. 



cit\0),c['^\0) 
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1) c^") 

^/^ ''22 



(1). Therefore the 
asymptotic factorization property is expected to be well- 
satisfied in this case. Moreover, the very weak cross- 
correlation between two systems in contact explains why 
the density profiles remain almost homogeneous even 
around the contact. Note that the density-correlation 



Jy) 



functions Ci^ (ri,^) and c{[ ifi,y) for the driven system 
along X and y directions, respectively, are clearly different 
due to the presence of a strong driving field which breaks 
the isotropy. Whereas, for the equilibrium system, cor- 
relations C22 (j'l.a;) and C22 {i"i.y) expectedly remain the 
same. For attractive interaction strengths, the correla- 
tions perpendicular to the driving field become negative 
at larger distances. This picture remains qualitatively 
the same at other densities as well. 



C. Excess Chemical Potential in the Driven KLS 

Model 

The role of the contact dynamics in the KLS model can 
be more complex than that in the previously discussed 
models of the ZRP. However, interestingly for some pa- 
rameter values, we indeed find evidence of a constant 
difference (5/i = /ii — /i2 in the excess chemical poten- 
tials in a wide density range even in the KLS model for 
nonzero driving. This supports the modified form of the 
large deviation principle even for the KLS model. 

We consider density n as a function of chemical poten- 
tial n for a driven system with K ~ —1, E ^ 6 which is 
separately kept in contact with various equilibrium sys- 
tems with interaction strengths K = 0, K = —0.8, and 
K = —0.9. As we have already discussed in Fig. [TTjfor 
the same systems, there were indeed deviations from the 
zeroth law observed as the different curves for n vs. fi 
in Fig. [TT] do not fall on each other. However, now by 
shifting chemical potential fi to jjl + 5^ and choosing a 
suitable Sfji in each of the cases, all the curves can be 
made to collapse on each other quite well as can be seen 
in Fig. [13 

We have done the same analysis for another driven 
system with a different set of parameter values K = Z.7b, 
E = 6. The driven system is separately kept in contact 
with two equilibrium systems with interaction strengths 
K = 0.75 and K = 1.0. In Fig. [161 the densities n are 
plotted as a function of the shifted chemical potential fi 
(also see Fig. [T^ . The curves could be collapsed on top 
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FIG. 15. Density n vs. shifted chemical potential /i is plotted 
for a driven system with K = —1, E = 6 which is separately 
kept in contact with three equilibrium systems with interac- 
tion strengths K = 0, K = -0.8, and K = -0.9. By shifting 
fj, ^ ^ + 5fi, all curves collapse on each other well where 
we choose 5fi — 0.06 and 5^ — 0.095 for the cases when the 
driven system is in contact with the equilibrium systems with 
K = —0.8 and K = —0.9, respectively. 




FIG. 16. Density n vs. shifted chemical potential fj. is plotted 
for a driven system with K — 3.75, E = 6 which is separately 
kept in contact with two equilibrium systems with K = 0.75 
and K = 1. Shifting fj, ^ fj, + S^, two curves collapse on each 
other reasonably well by choosing Sfj, = 0.1 for the case when 
the driven system is in contact with the equilibrium system 
with K = 0.75. 



of each other reasonably well by shifting the fi to fi + 5fi. 
However, in both the cases of Figs. [TS] and [TBI it should 
be noted that the collapse is not so good at very low 
chemical potentials (i.e., low densities). This indicates 
that the difference in the excess chemical potential dfi 
actually may not be constant over the entire chemical 



potential range and can depend on the chemical poten- 
tial itself (or equivalently density) of the corresponding 
equilibrium system in contact. 



VI. CONCLUDING PERSPECTIVE 

In summary, in this paper we have studied "equilibra- 
tion" of driven lattice gases when two systems are kept in 
contact and allowed to exchange particles with the total 
number of particles conserved. Interestingly, both for at- 
tractive as well as for repulsive nearest-neighbor interac- 
tions and for a wide range of parameter values, there is a 
remarkably simple, though approximate, thermodynamic 
structure where the zeroth law is quite well satisfied. 

However, there are also observed deviations from this 
simple thermodynamic law. To understand these devi- 
ations, first we have studied the nontrivial role of the 
contact dynamics for a variant of the zero range pro- 
cess (ZRP) as well as a variant of the equilibrium Katz- 
Lebowitz-Spohn (KLS) model, where one can calculate 
the steady-state probability distribution exactly. Using 
these simple examples, we point out that, due to the 
modified contact dynamics, there can be an excess chem- 
ical potential induced across the contact region. These 
results lead us to express the large deviation principle in 
a modified form which elucidates the role of the contact 
and which can substantially account for the deviations 
from the zeroth law. 

It is important to note that the modified form of the 
large deviation principle is still based on the asymptoti- 
cally factorized form of the steady-state distribution for 
two driven systems in contact and is valid if the cor- 
relations between the systems can be neglected in the 
large volume limit. In the case of the variant of the ZRP 
discussed here, these correlations are zero. In the case 
of the driven KLS model, we observed that the spatial 
density-correlations across the contact are indeed very 
small compared to the nearest-neighbor correlations in 
the individual systems. Interestingly, we found evidence 
of an almost constant excess chemical potential for var- 
ious parameter values and in a wide range of densities, 
therefore supporting applicability of the modified large 
deviation principle to the KLS model. 

In general, the results in this paper lead us roughly 
to the following possible scenarios for driven systems in 
contact. 

(1) The large deviation principle, or in other words the 
asymptotic factorization property, may break down due 
to long-ranged correlations which may be present in the 
driven systems |2£ 



27[. In this case, the combined sys- 



tem cannot be divided into independent subsystems and 

there would be no intensive variable which equalizes upon 

contact. 

(2) However, when the amplitude of these correlations 

are sufficiently weak, it is possible that a large deviation 
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principle holds, although in a modified form. In these 
cases, the systems can be characterized by the excess 
chemical potentials across the contact and consequently 
there would be some intensive variable which would then 
equalize upon contact. Note that introducing the excess 
chemical potential is essentially a way of reparametriza- 
tion of the chemical potential of the driven system under 
consideration. These excess chemical potentials can de- 
pend on the specifics of the contact dynamics and they 
are generally a-priori unknown, therefore in a sense ar- 
bitrary. For some parameter values, the excess chemical 
potential may almost be constant over a range of den- 
sities. In special cases, the arbitrariness of the excess 
chemical potential can be removed by choosing a suit- 
able contact dynamics such that the zeroth law holds 
strictly. In these cases, it is actually possible to assign 
to the individual systems an intensive variable, indepen- 
dent of the specifics at the contact between two systems, 
which would then equalize. It should be noted here that 
the modified large deviation principle could be satisfied 
irrespective of whether the zeroth law is satisfied or not, 
which has been illustrated in this paper by using a variant 
of the ZRP. 

For the ZRP, since the steady state properties are ex- 
actly known, it is easier to choose a contact dynamics so 
that the zeroth law can be made to be satisfied. How- 
ever, for the KLS model which has nontrivial steady state 
properties, it would be difficult to find a contact dynam- 
ics, even if it exists, for which the zeroth law would hold 
strictly. For the KLS model, here we have considered 
a contact dynamics, mainly based on a physical ground 
albeit on an ad-hoc basis, which satisfies the local de- 
tailed balance condition. The simple modification of the 
large deviation principle suggested in this paper indicates 
that it could still be possible to choose a contact dynam- 
ics such that the thermodynamic laws are satisfied even 
better. As an open question, the origin of the excess 
chemical potential at the contact should be understood 
in more detail which could give valuable insights whether 
it is possible to choose a contact dynamics for which a 
simple thermodynamic structure emerges for driven sys- 
tems in general. 

There are further important aspects in exploring such 
a simple structure for driven systems in contact. Unlike 
in equilibrium where there is a well defined prescription 
to describe various thermodynamic properties using the 
standard Boltzmann distribution, there is no such pre- 
scription for nonequilibrium systems. However, the nu- 
merically observed simple thermodynamic structure con- 
cerning driven lattice gases in contact may give us an 
useful tool to characterize such systems and gives rise 
to a possibility to describe phase transitions which are 
known to occur in these driven interacting many-particle 
systems. 
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APPENDIX A: STEADY STATE DISTRIBUTION 
FOR THE ZRP WHEN ui / V2 



In the general case when vi 



„W 



^ V2 



^2 > 



the ansatz for the steady-state probability distribution 
is given by 



Pst{{n^^},{n^^}) = -^ 



\{h{n^,)^h{^ 



1 



iVi ,JV2 



V-i V. 



12 = 1 

■6{Ni +N2-N) (19) 



where A^i and N2 are number of particles in ring 1 and 
ring 2 respectively and A^ = A^i -I- A^2 is the total con- 
served particle number. Now we consider the following 
two cases - (1) when particle jumps in the bulk, say in 
ring 1, and (2) when particle jumps from one ring to the 
other. 

Case (1): Consider the following two transitions 
from a configuration C to C" and a configuration 
C" to C where C = {{. . . ,ni^-i,ni^,n.i^+i, . . ■},{ni^}), 
C = {{. . . ,ni^^i,n.i^-l,n.i^+i + l, . . .},{ni^}) and C" = 
({..., ni^_i -I- l,ni^ - l,riij+i,...},{rij2}). The steady- 
state probability current J{C' -^ C") = Pst(C)w(C"|C) 
from configuration C to C", Pst{C) being the probability 
of configuration C in steady state and w{C'\C) being the 
transition rate from C to C", can be explicitly written as 

J{C ^ C) = ^ [. . . /K_i)/(n,J/(n,,+i) . . .] 



(b)fint, - 1) 



1 

■ Ni N2 



11 /2KJ 



.i2 = l 

5{Ni +N2~N) (20) 



where we have used jump rate in the bulk of ring 1 which 
is w{C'\C) = vf^f{n,, - l)//(n,J. Finally we get 

■AC ^ C') = ^ [. . . /(n.,-i)/(n,, - l)/(n.,+i) . . .] 



Vih{^ 



.i2 = l 



V 



(6) 



A^i„,JV2 



Wi V. 



S{Ni+N2-N). (21) 



Similarly the probability current J{C" -^ C) ~ 
P{C")w{C\C") from configuration C" to C can be ex- 
plicitly written as 



J{C" ^ C) - — [. . . /(n,,_i + l)/(n,, - l)/(n,,+i) . 
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fir 



'ti-D 



(b) 

"' /K-i + 1) 
1 



Uhi^ 



S{Ni +N2~N) 



Zn 

' L2 

n /2K) 



/("-Jl-l)/("-ll - l)/("-Jl + l) ■ • •] 



.ib) 



..Ni„.N: 



■S{Ni+N2-N). (22) 



Clearly the probability currents J{C — > C") and 
J{C" -^ C) are equal. Since, for any transition from 
C to C, it is always possible to find a corresponding 
transition C" to C, the net current into C vanishes 
pairwise, i.e., J{C -^ C) ~ J{C" -^ C) in the steady 
state. 

Case (2): Consider two transitions from C 

to C" and the configuration C" to C where 
C = ({...,nfci,...},{...,nfe2, ...}), C" = 
({..., Tifc^ - !,...},{. ..,nk2 + 1, ■•■}), and fci and 
/c2 are respective contact sites in ring 1 and ring 2. In 
this case the probability current J{C — > C") can be 
written as 



J{C -^ C") 



^[■■•/KJ---] 



Vl- 



X [. ../(nfcj.. 






1 

Zn 



■fink, -I). 
1 



Wi-l ^2 



w 



Vn 



/Kl) 

<5(iVi + 7V2 - TV) 



(5(iVi + 7V2 - iV)(23) 



Similarly the probability current JiC' — !• C) can be writ- 
ten as 



J{C' 






./K, - !)...][... /K,+l) 



W2 



/K2 + 1) 



1 



,^1 



-1 N.+A ^i+N2-N) 



1 



1 



1[.../KJ...] 

" Ar,_i ^/(iVl + A^2 - iV). (24) 
Vi f 2 

Clearly the net probability current into C is again zero 
as J{C -^ C) = J{C' -^ C) in the steady state. This 
completes the proof for the steady-state ansatz given in 
Eq. [T^ which satisfies Master equation dtPstiC,t) = = 
Ec'^c[PstiC')wiC' ^ C) - PstiC)w{C -^ C')]. Note 
that the condition J{C -^ C) = J{C' -^ C) is nothing 
but the detailed balance condition which is satisfied at 
the contact even when vi ^ vi- However, detailed bal- 
ance is not satisfied in the bulk and consequently there 
are nonzero currents within the individual systems. 



APPENDIX B: FLUCTUATION-RESPONSE 
RELATIONS 

One interesting consequence of Eq. [15] is a relation be- 
tween the susceptibility and the fluctuation in particle- 
number of a system when the system is in contact with 
a large reservoir characterized by a chemical potential /x. 
Consider system 1 is in contact with system 2 which is 
very large compared to system 1. Let us denote a^^^ as 
the standard deviation of fluctuations in the total num- 
ber of particles A^i in system 1, i.e., cr^ = (iVf) — (iVi)^. 
Then the large deviation principle with the definition 
of chemical potential as given in Eq. [15] implies that 
the fluctuation of particle-number around an average 
particle- number A^*, 



P(7Vi) 



constant x e 



-{N,-Nlfl2x 



(25) 



where 1/x = d^ildN^ = {d^Fi/dN^N^ with Fi{Ni) = 
Vifi{ni). Since the root mean square fluctuation in the 
particle-number A^i is cr^ = 2x, one gets the following 
fluctuation relation, 






= ((iV^) - {N,f). 



(26) 
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FIG. 17. Driven system in contact with an equilibrium reser- 
voir of noninteracting hardcore particles. Integrated com- 
pressibility I-), (squares) and integrated fluctuation I^ (circles) 
plotted as a function of fi for a driven system with two diiTer- 
ent sets of parameter values with K = —1, E = 6 (red) and 
K = —0.75, E — 4 (blue). The fluctuation-response relation 
is well satisfled. 



For the ZRP, even when vi ^ V2, the fluctuation- 
response relation between the compressibility and the 
root mean square fluctuation in particle-number A'^i is 
still exactly satisfied if the ring 1 is in contact with a 
much larger ring 2 being a particle reservoir. Moreover, 
since in this case the variables 11'^ and fj^a differ only by 
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a constant In Va (see Eq. [10] and the paragraph below) , 
the fluctuation-response relation is satisfied with respect 
to both the new and old intensive variables n'^ and Ha, 
respectively. 

Now wc briefly discuss the numerical results concern- 
ing the fluctuation-response relation for the KLS model 
with attractive interaction strengths. To numerically test 
the fluctuation relation as given in Eq. 1261 we consider 
a driven system which is in contact with an equilibrium 
reservoir of hardcore particles which are otherwise non- 
interacting (i.e., K = Q). The chemical potential /x of 
the equilibrium hardcore particle-reservoir is given by the 
expression in Eq. |4l For better numerical accuracy, we 
check the integrated version of Eq. [26l i-e., we calculate 
the integrated fluctuation /„ = /^ (cr^Jd/i' and the inte- 
grated susceptibility /^ = /'' {d{Ni) / d^')d^' for different 
values of /i obtained by varying density of the equilibrium 
reservoir. We take a two-dimensional 20 x 20 nonequilib- 
rium system in contact with a 250 x 250 equilibrium reser- 
voir. In Fig. I17[ the integrated compressibility and the 
integrated fluctuations are plotted as a function of chem- 
ical potential /i for two driven systems with two different 
sets of parameter values with [1) K = —\, E = Q and 
(2) K = —0.75, E = A. In these cases, the fluctuation- 
response relation is remarkably well satisfied as seen in 

Fig. [m 
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